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I investigate the phenomenology of supersymmetric models with extra vector- 
like supermultiplets that couple to the Standard Model gauge fields and trans- 
form as the fundamental representation of a new confining non-Abelian gauge 
interaction. If perturbative gauge coupling unification is to be maintained, the 
new group can be SU(2), SU{3), or 5*0(3). The impact on the sparticle mass 
spectrum is explored, with particular attention to the gaugino mass dominated 
limit in which the supersymmetric fiavor problem is naturally solved. The new 
confinement length scale is astronomical for 5*0(3), so the new particles are 
essentially free. For the SU (2) and SU (3) cases, the new vector-like fermions 
are quirks; pair production at colliders yields quirk-antiquirk states bound by 
stable fiux tubes that are microscopic but long compared to the new confine- 
ment scale. 1 study the reach of the Tevatron and LHC for the optimistic case 
that in a significant fraction of events the quirk-antiquirk bound state will lose 
most of its energy before annihilating as quirkonium. 
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I. INTRODUCTION 

Among the hurdles that must be cleared by any proposed extension of the Standard 
Model (SM) are the stringent limits on quantum corrections to the electroweak vector boson 
propagators due to new physics [H-jsl. Low-energy supersymmetry [3] is generally safe in 
this regard, because of the fact that all of the new particles it introduces get their masses 
primarily from bare mass terms, not from their couplings to the Higgs vacuum expectation 
values (VEVs). This includes the Higgs chiral supermultiplets and themselves, which 
are vector-like, together forming a self-conjugate representation of the SM gauge group 
GsM = 5't/(3)c X SU{2)l X U{1)y- It is therefore interesting to consider non-minimal 
supersymmetric models that maintain this feature by including more chiral supermultiplets 
transforming as vector-like representations of the gauge group. 

Another well-known and appealing feature of the minimal supersymmetric standard 
model (MSSM) is the perturbative unification of running gauge couplings near Mu 2x 10^^ 
GeV. A sufficient (but not necessary) condition for extensions of the MSSM with extra 
vector-like supermultiplets to maintain gauge coupling unification is that the new fields 
come in complete multiplets of the SU{b) global symmetry group that contains Gsm- This 
paper studies the properties of models of this type that introduce a new non-Abelian gauge 
group Gxi under which the new chiral supermultiplets also transform but the MSSM fields 
are neutral. Models of this type have already been introduced by Babu, Gogoladze, and 
Kolda in jsl, in the context of finding new large contributions to the lightest Higgs boson 
mass.^ I will assume that the new chiral supermultiplets have masses at the TeV scale or 
below, and that the new gauge coupling gx unifies with the Gsm couplings gi, g2, gs at Mu. 
In order to avoid a strong disruption of the running of the Gsm gauge couplings, it is neces- 
sary that the corresponding confinement scale A for the Gx interactions is below the masses 
of the new fermions and scalars that are also charged under Gsm- This in turn implies an 



intriguing phenomenology studied first by Okun [15|, later by Gupta and Quinn [16| and 



by Strassler and Zurek [17[, and more recently in considerable depth by Kang and Luty in 
18|. The new particles that transform non-trivially under the new gauge group (dubbed 
"theta particle" by Okun, and renamed "quirks" by Kang and Luty) can form exotic bound 
states with unusual signatures that depend strongly on the Gx confinement scale. When a 
heavy quirk- ant iquirk pair is produced in a collider experiment, they fly apart but remain 
connected by a stable flux tube, which cannot break due to the large energy cost to produce 
an additional quirk- ant iquirk pair. The maximum length of this flux tube is roughly of or- 
der L ~ AE/A"^, where AE is the kinetic energy of the hard scattering production process. 
This length can range from microscopic to literally astronomical, but in any case it is much 
larger than the flux tube thickness ~ A~^. The resulting collider signatures are potentially 



distinctive but also possibly quite difficult [l8|-[27 . 

In this paper, I will study the basic properties of models that maintain perturbative 
unification of gauge couplings, and their renormalization group running, in Section [Tll The 



sparticle mass spectra are studied in Section [Hit and Section [TVl considers the impact on the 



t Models with the same motivation, but without the new non-Abehan gauge group, have been studied in 



13| . Other recent proposals for extra vector- like chiral supermultiplets are found in 14 
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supersymmetric little hierarchy problem. Some salient aspects of the collider phenomenology 
of the quirks are discussed in Section |Vl 

II. MSSM EXTENDED BY VECTOR-LIKE FIELDS COUPLED TO A NEW 
CONFINING NON-ABELIAN GAUGE INTERACTION 

In order to maintain perturbative gauge coupling unification, the number of new particles 
transforming under the SM gauge group is limited to the equivalent of three copies of the 
5 + 5 of the SU{5) group that contains Gsm, if they are not much heavier than 1 TeV. This 
assumes that = gf /in {i = 1, 2, 3) are required to be perturbative (less than 0.3 or so) at 
and below the energy scale where they unify.-f (One could consider unification with larger 
couplings at and near the unification scale, but then both renormalization group (RG) 
running and threshold corrections will be necessarily out of control, and the low-energy 
manifestation of apparent unification must be considered merely accidental.) It follows 
that the new gauge non-Abelian group must be Gx = SU{2)x or S0{3)x or SU{3)x- In 
the following, the new fields are taken to transform in the = 2, 3, or 3 dimensional 
representations respectively for these three cases. Thus the new quirk chiral supermultiplets 
transform under SU{2)x x Gsm as: 



D,D = (2, 3,1,-^) + (2, 3,1,^; 

L,L = (2,l,2,i) + (2,l,2,-^; 
S,S = (2,l,l,0)x2n5, 

or under SU{3)x x Gsm as: 

D,D = (3, 3,1,-^) + (3, 3,1,^) 

X,L = (3, 1,2,^) + (3, 1,2, 



2' ' ' ' 2' 
S,S = [(3,1,1,0) + (3,1,1,0)] xn5. 



or under SO{?>)x x Gsm as: 



D,D = (3, 3,1,-^) + (3, 3,1,^) 

L,L = (3,l,2,i) + (3,l,2,-i) 
S = S = (3,1,1,0) X ns. 



2.1 

2.2 
2.3 

2.4 

2.5 
2.6 

2.7 

2.8 

2.9 



These are the main model frameworks considered below. With these assignments, L 
transform as a 5 and Z), L transform as a 5 of the usual Georgi-Glashow SU{b), ensuring 



^ It is crucial to use two-loop (or higher) beta functions to correctly implement this perturbativity require- 
ment. This paper uses three- loop beta functions for supersymmetric gauge couplings and gaugino masses 
and two-loop beta functions for Yukawa couplings, scalar masses, and scalar cubic couplings. These can 



be found straightforwardly from general results in refs. j28l |29|. and so arc not listed explicitly here. 
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that the unification of Gsu gauge couphngs persists. Also included are ns Gsm singlets in 
the same representations of Gx- Note that since SU{2)x and SO{?>)x have the same Lie 
algebra, the practical distinction between them is really whether the representations of the 
chiral superfields are doublets or triplets. 

[There are some variations on the above models that are consistent with gauge coupling 
unification with the new fields at the TeV scale, which should be mentioned although they 
are inconsistent with an assignment of -D, -D, L, L into 5 + 5 of SU{5). First, for SU (3)x only, 
there is another, inequivalent, embedding in which D, D have the same assignments, but 
L,L = (3, 1, 2, i) + (3, 1, 2, — i) instead. Also, for SU{3)x, one could put D, D into three 
singlets if L, L are in a triplet, or vice versa. Likewise, for SU{2)x or S'0(3)x, one could put 
D, D into any combination of nf singlets, doublets, and triplets, and similarly for 
L, L, provided that nf + 2nf + 3^1^ = nf + + 3^3 < 3. Finally, it should be noted that 
the number and type of Gx representations of the SM singlets do not affect gauge coupling 
unification for Gsm, and so are more generally arbitrary as long as they are anomaly-free 
under Gx- However, to keep the discussion below bounded, I will limit the discussion below 
to the models defined by eqs. ([2l])-([23D.] 

The supersymmetric mass parameters of the D, L, L fields are assumed to arise by the 
same mechanism that gives the entirely analogous term fiHuHd in the superpotential of the 



MSSM. For example [30|, l3l|, one may assume that the mass terms HuHd and DD and LL 



and SS are forbidden at tree-level, and arise from non-renormalizable superpotential terms: 
W = -^XX + XnDD + XlLL + A5, 5.5,) (2.10) 

(with an implied sum over 2 = l,...,n5ifn5 7^0) when the fields X, X get VEVs roughly of 
order 10^^ GeV. Here Mp = 2.4 x 10^^ GeV is the reduced Planck mass. These intermediate- 



scale VEVs are natural, for example 3l|, if there is also a superpotential 



W^^^X^X (2,11) 



and soft terms 



- Aoft = m\\X\' + m||X|2 + i^-^X'^X + c.c. j . (2.12) 

Non-trivial VEVs for X, X break a Peccei-Quinn symmetry, giving rise to an invisible axion 
solution to the strong CP problem [3^. There will be a non-trivial local minimum of the 
potential provided that |axP — 6 1 A y P (my + m^) > 0, and it will be a global minimum if 



o-xP — 8|AxP('^x + ^ ^ This will give rise to the vector-like mass terms in the 



low-energy effective superpotential 

W = fiHuHd + fioDD + fiLLL + fisASi- (2.13) 
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with jj,, Ho-, jJiL^ l^s of order 100 GeV to 1 TeV, provided that the corresponding couphngs 
A/i, A/), Ai, \s are not too small 

The S, S fields do not couple to Gsm gauge fields, and so are not constrained by LEP2 
or Tevatron or other direct production, nor do they affect the SM gauge couplings directly. 
So, some number n of them (with < n < ns) could actually have current masses ixsi that 
are far below the electroweak scale. This would occur if the A^^ coupling(s) in eq. f l2.10p are 
absent (perhaps replaced by terms of even higher dimensionality), or just small. Note that 
for n > 0, there will be no stable flux tubes for pair-produced particles charged under Gx-, 
because then as the particles produced in the hard collision fly apart, the gauge string will 
break to form bound states with size of order just as in ordinary QCD. This is because 
the energy cost to produce an additional pair of light S*, S after the hard collision would 
then be small. 

The new fermion content of the theory consists of a color triplet charge ±1/3 Dirac 
fermion r, , ifj j-,) with mass /xu; a charge ±1 Dirac fermion {ip^,ip2) with mass /i^; and 
charge fermions . The scalar partners of these particles will have soft- 

supersymmetry breaking squared-mass terms: 

- £ = m||D|2 + m||D|2 + m||L|2 + ml\L\'^ + {ml)ijS*Sj + {m'^)ijS*Sj + 

+ {bDDD + bLLL + {bshjSiS, + c.c), (2.14) 

where the scalar components are denoted by the same symbol as the chiral supermultiplets 
of which they are members. In the case Gx = SU{2) or 5*0(3), the flelds Si and Si actually 
have the same quantum numbers, and so can mix with further soft mass terms {■m~)ijS* Sj, 
etc., but for simplicity I assume that mixing between Si and Sj chiral supermultiplets is 
absent. Also for simplicity, I will assume that the above Si and Sj soft terms are diagonal 
in the same basis that the superpotential masses are diagonal. This is natural if the soft 
supersymmetry breaking arises in a flavor-blind framework such as gaugino mass dominance. 

In the absence of Yukawa couplings involving the new chiral supermultiplets, the charge 
fermions are unmixed, and form Dirac fermions with masses hl and fis,- For ns > 0, the 
new chiral supermultiplets can have Yukawa couplings in addition to their mass terms in 
eq. (EISD: 

W = kiH^LSi + KHdLS,. (2.15) 

and corresponding soft scalar cubic terms, 

- C = ak^HuLSi + ay.HdLSi + c.c. (2.16) 

As mentioned above, H^LS and H^LS couplings for SU{2)x or S0{3)x are also possible, 
but are omitted here for simplicity. The superpotential Yukawa couplings produce mixing 
between the gauge eigenstate fermions, yielding Dirac fermions {ip^,ipj) which are mixtures 
of ipl, ijjg^ and of ip^, , respectively. For example, if only one pair 5* and S has couplings 
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to the Higgs fields, then the new neutral fermion mass matrix becomes 



-C = I ^^"1 I - I +C.C. (2.17) 




when the MSSM Higgs fields get their VEVs f^, with tan/3 = v^/vd and v = v^t^^Tuf ~ 
175 GeV. The couplings k, k' gives rise to 1-loop effects that can significantly raise the 
lightest Higgs scalar boson mass due to a lack of complete cancellation between scalar and 
fermion loops, especially for large k if tan/3 is not small. This was the motivation of jsf. 



but as noted in similar contexts in |10l . Illj and remarked on further below, it is doubtful 
whether this really ameliorates the supersymmetric little hierarchy problem. 

In keeping with the idea that the apparent gauge coupling unification for Gsm is telling us 
something important about the underlying theory, I will assume that the new non-Abelian 
gauge coupling gx unifies with gi, g2 and g^ at a scale Mu > 2 x 10^^ GeV. In practice, 
I use three-loop RG equations to run up from the electroweak scale, and declare the scale 
where gi = g'^y5/3 and g2 meet to be Mu, and require gx to be equal to them there. The 
QCD coupling g^ typically misses this common value at Mu by a small amount that can 
be reasonably ascribed to threshold corrections. Now, RG running gx from this scale, I 
require that it remains finite down to scales well below the masses of the quirks D, D, L, L. 
Otherwise, two-loop effects would strongly affect the running of the SM gauge couplings, 
rendering their apparent unification merely accidental. For SU{2)x and S0{3)x, this re- 
quirement is automatically satisfied for all ns > 0, but for SU{3)x it requires ns > 3. 1 
therefore consider ra^ = 3 to be the minimal viable model for the SU {3)x case. 

For illustration, the running of the gauge couplings is shown at three-loop order in Figure 
1, for the three cases SU{2)x with ns = and S0{3)x with ns = and SU{3)x with 
ns = 3. For simplicity, I have assumed vanishing Yukawa couplings and chosen a single scale 
M;hresh = 1 TeV as the effective average mass of the new particles charged under Gx and the 
MSSM superpartners. The unification will have some dependence on the actual thresholds, 
which one might imagine is roughly comparable to the unknown threshold dependence due 
to high-scale particles. In the case of SU{2)x with ns = 0, the gauge coupling gx runs quite 
slowly, and is somewhat weaker than the QCD coupling at the TeV scale. In contrast, for 
the S0{3)x case with ns = 0, gx runs quickly to very small values in the infrared, due to 
a large positive beta function coefficient. For the minimal viable SU{3)x case with ns = 3, 
the gx beta function is even more negative than the QCD beta function, leading to a gauge 
coupling at the TeV scale that is larger, but still perturbative and not running very fast. 
For non-minimal models with ns larger than these values, the TeV-scale values of ax are 
smaller, because the gx beta function is larger. 

Below the masses of the quirks and their supersymmetric partners, the coupling gx has 
a negative beta function, and diverges at some scale A when calculated at any particular 
loop order in a specified scheme. Given the MS beta function for ax = 5'l:/4vr up to 4-loop 
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FIG. 1: Renormalization group running of the 
inverses of the gauge couphngs ai, 02, 03, and 
ax, for the minimal viable cases Gx = SU{2) 
with ns = (upper left), Gx = SU{3) with 
71-5 = 3 (upper right) and Gx = SO (3) with 
ns = (lower left). The undiscovered MSSM 
particles and the new particles charged under 
Gx are taken to have masses at a common 
threshold Mthresh = 1 TeV. 



order: 



Pa, = Q- 



da 



X 



dQ 



-2 (6oax + + hax + ha% + . . .) 



(2.18) 



the scale A can be defined,^ using any convenient ax{Qo) with Qo < Mthrcsh as input, by an 
expansion in inverse powers oi t = \n{Ql/A'^) [33| : 

«x(go) = ^(l-[h\nt]/blt+[boh + bl{lnH-\nt-l)]/btt^ 

OqI \ 

+ [6^63 - &?(2 In^ t - 5 In^ t - 4 Int + 1) - In t\/2blt^ + ...). (2.19) 

It is common in rough estimates to only use the one-loop-order estimate ax{Qo) = ^/bot, 
with bo = {UCa - 2rp)/127r where (Ca,Tp) = {N,n) for SU{N)x and {Ca,Tf) = (2, 2n) 
for S0{3)x, with n denoting the number of the ns SM singlet fields that have masses below 



§ Note that the definition for A used here corresponds to A/4 in ref. 



181 
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Gx 


ns 


TeV) 


Ad) 


A(2) 


A(3) 


A(4) 


SU{2) 





9.3 


0.35 GeV 


1.3 GeV 


1.1 GeV 


1.1 GeV 




1 


14.4 


4.4 MeV 


19 MeV 


17 MeV 


17 MeV 




2 


19.5 


57 keV 


280 keV 


250 keV 


250 keV 




3 


24.5 


0.76 keV 


4.1 keV 


3.7 keV 


3.7 keV 




4 


29.5 


11 eV 


60 eV 


55 eV 


55 eV 


SU{?,) 


3 


4.9 


61 GeV 


140 GeV 


120 GeV 


120 GeV 




4 


9.9 


3.5 GeV 


11 GeV 


9.3 GeV 


9.5 GeV 




5 


15.0 


190 MeV 


720 MeV 


620 MeV 


620 MeV 




6 


20.1 


10 MeV 


44 MeV 


38 MeV 


38 MeV 




7 


25.1 


0.59 MeV 


2.7 MeV 


2.4 MeV 


2.4 MeV 




8 


30.2 


32 keV 


160 keV 


140 keV 


140 keV 




9 


35.2 


1.9 keV 


9.7 keV 


8.7 keV 


8.8 keV 


50(3) 





83 


1.3x10^19 eV 


1.1x10-1^ eV 


l.OxlO-is eV 


1.0x10-18 eV 




1 


103 


4.7x10-27 eV 


4.4x10-26 eV 


4.2x10-26 eV 


4.2x10-26 eV 



TABLE I: The Gx confinement scale A^^, computed at various loop orders I by using eq. (|2.19|) 
keeping terms of order \/t^. The new particles charged under Gx are taken to have an effective 
average decoupling mass scale of Qq = Mthresh = 1 TeV, with no SM singlet quirks with current 
masses smaller than A. 

Qo and are treated as non-decoupled below Mthresh- However, it turns out that including 
the higher loop effects (with coefficients fei,2,3 found in refs. (34]) are quite important for 
obtaining a stable value of the Gx confinement scale A. This is illustrated in Tabled which 
shows the results obtained for A^ at various loop orders assuming again that the effective 
decoupling scale for particles charged under Gx is Mthresh = 1 TeV. The point of carrying 
the calculation to 4-loop order is not because of the very slightly increased accuracy obtained 
(since there are threshold uncertainties here that are not known), but rather to demonstrate 
the stability of the results with respect to inclusion of higher-order terms. In fact, the 4-loop 
order results for A hardly differ at all from the 3-loop order ones, and only at the 10% level 
from the 2-loop order ones. However, they are notably larger than the 1-loop order estimate, 
which is therefore judged to be deprecated as an estimate of the physical Gx confinement 
scale. 

Table [T] shows that the confinement scale A for S'0(3)x is very small in energy units. In 
terms of length, the confinement scale for the minimal model = is of the order 10^^ 
meters, very roughly of order the radius of the Earth's orbit around the Sun. For = 1, 
the confinement length is of order 100 parsecs. Thus for all practical purposes, the quirks 
are actually free. Adding other SM singlets charged under SO{?>)x will only decrease A, 
making the confinement length even larger. 

For the minimal viable SU{2)x and SU{3)x models, the confinement energy scale is 
much larger. Increasing ns leads to smaller A, as indicated in Table [H If n of the ns SM 
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Gx 


ns = n 


ax\l TeV) 


A(4) 


SU{2) 


1 


14.4 


5.0 MeV 




2 


19.5 


6.3 keV 




3 


24.5 


1.3 eV 




4 


29.5 


2.0 X 10-5 eV 


SU{3) 


3 


4.9 


68 GeV 




4 


9.9 


1.5 GeV 




5 


15.0 


13 MeV 




6 


20.1 


38 keV 




7 


25.1 


30 eV 




8 


30.2 


0.0030 eV 




9 


35.2 


2.1 X 10-s eV 



TABLE II: As in Table [H but now taking all ns = n oi the SM singlet fermions to have current 
masses much less than A, and showing only the 4-loop order result (which is nearly identical to 
the 3-loop order result in all cases). 



singlets charged under Gx have current masses /i^. less than A, the confinement scale will 
be decreased. This is illustrated in Table |TT] for the extreme case that all n = ris oi the 
new singlets are lighter than A. Both Tables HI and HTl take the effective average decoupling 
scale for the particles in the other new chiral supermultiplets (including both scalars and 
fermions) to be 1 TeV. More generally one can estimate: 

/ ^thresh A ^ ^ . 3Cyd — 5 — nc , , 

A = A..,^^j ■ ^- (llC.-2n)/3 - '^■^»> 

where now Atab is the value given in Table |T] for n = 0, or Table HTl for n = ns- Here, Mthresh 
is defined to be the effective average decoupling scale for the new supermultiplets, at which 
the threshold corrections to the gauge coupling gx are small. As seen in Figure [H in the 
minimal viable models the gx coupling runs fairly slowly in the non-decoupled theory above 
Q = Mthresh- This mcaus that in the minimal model for SU{2)x, A ^ O.OOlMthresh, while 
for the minimal viable SU{3)x model with = 3, A 0.12Mthresh, if n = 0. 

If the Yukawa couplings k, k' are present, there is a potentially important constraint 
from precision electroweak observables. The new contributions to the Peskin-Takeuchi S, T 
observables from the new fermions are: 

F 
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where k = and k' = fc'cos/3 and f ^ 175 GeV, and for illustration purposes I have 

chosen fiL ^ f^s ^ and assumed that the corresponding scalars are much heavier. The 
values of these Yukawa couplings are governed by infrared quasi-fixed points. For example, 
if k' is negligible, then the beta functions for k and the top Yukawa coupling are given at 
one-loop order by: 



Q— 
^dQ 



(3k 



(3 + + 3y^ - ACgl - ?,gl - ^1 



(2.23) 
(2.24) 



where iV = 2, 3, or 3 and C = 3/4, 4/3, or 2 for Gx = SU{2), SU{3), or ^0(3) respectively. 
The fixed points arise due to the balancing between the positive Yukawa and the negative 



gauge contributions 35|. Including two-loop effects, I find for the minimal viable models 



the infrared quasi-fixed-point values: 



''fixed 




[SUi2)x,ns = 0] 
[SO{3)x,ns = 0] 
[SU{3)x,ns = 3]. 



(2.25) 



at Q = 1 TeV. The resulting contributions to S, T can be used to put a lower bound on 
fiL- Requiring the results to be within the current 95% CL ellipse from experimental results 



on rrit, mw, and Z-peak observables using the same methodology as in I estimate 
fiL > 210,225,380 GeV for the Gx = SU{2), S0{3), SU{3) fixed point cases respectively. 
However, Gx confinement may play a significant role in modifying this estimate for SU{3), 
because A in that case is larger than Mz- For smaller Yukawa couplings k <^ kfi-^^d, there is 
no constraint as the vector-like particles decouple from precision electroweak observables. 



III. SOFT SUSY-BREAKING MASSES AND THE SPARTICLE SPECTRUM 

The presence of new vector-like supermultiplets has a profound effect on the spectrum of 
superpartner masses. They cause the Gsm gauge couplings to run to much larger values in 
the ultraviolet as they approach unification, resulting in bigger one-loop contributions to soft 
scalar squared masses from RG running, compared to the MSSM. The new supermultiplets 
also allow the Gx gaugino masses to contribute indirectly to MSSM gaugino and sfermion 
masses, through two-loop order effects. In this section, the patterns of soft supersymmetry 
breaking masses will be considered for these models. For simplicity, the discussion will be 
mostly limited to the scenario in which a unified gaugino mass parameter mi/2 is much larger 
than the scalar masses and other sources of supersymmetry breaking at the RG scale where 
the gauge couplings unify. This gaugino mass dominated limit is motivated as a solution to 
the supersymmetric flavor problem, since it automatically produces flavor-blind soft terms. 

The modified running of the gaugino masses pushes them to be smaller near the TeV 
scale than they would be in the MSSM. Given an input unified gaugino mass mi/2 at the 
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unification scale, one finds for the running gaugino masses at Q = 1 TeV: 



(Ml, M2, M3) = mi/2 X < 



(0.41, 0.77, 2.28) 
(0.21, 0.39, 1.16) 
(0.112, 0.20, 0.57) 
(0.080, 0.135, 0.40) 



[MSSM], 
[SU{2), ns = 0], 
[50(3), ns = 0], 
[SU{3),ns = 3]. 



(3.1) 



Thus, in the extended models, to obtain the same physical gaugino masses, one must start 
with larger mi/2 than one would in the MSSM. Since mi/2 is not directly observable, it is 
also interesting to consider the ratios of these gaugino masses. They are also affected, but 
more mildly (being due to 2- loop effects): 



(M2/M1, M3/M2, M3/M1 



^1.87, 2.96, 5.53) 

;i.85, 2.95, 5.44) 

;i.79, 2.84, 5.09) 

;i.73, 2.85, 4.95) 



[MSSM], 
[SU{2),ns = 0], 
[SO{3),ns = 0], 
[SU{3),ns = 3], 



(3.2) 



where again unification of gaugino masses at the gauge coupling unification scale is assumed. 
The effect of the additional fields is thus to somewhat compress the gaugino mass spectrum 
compared to the MSSM case, with the ratio of gluino to bino masses decreased by about 10 
per cent for Gx = 5*0(3) and SU{3). To obtain the physical masses, one must also include 
mixing with Higgsinos and the pole mass corrections, which are particularly important for 
the gluino 36l-l38l . 

In the extended models the squark and slepton masses are also relatively smaller (com- 
pared to mi/2) at the TeV scale than in the MSSM. Taking a gaugino-mass dominated 
scenario (by assuming a vanishing common scalar squared mass = at the unification 
scale), one finds for the first and second family squark and slepton masses at Q = 1 TeV: 



mi/2X < 



(2.15, 2.08, 2.07, 0.67, 0.37) 
(1.61, 1.55, 1.54, 0.57, 0.33) 
(1.23, 1.18, 1.17, 0.48, 0.28) 
(1.06, 1.02, 1.01, 0.43, 0.26) 



[MSSM], 
[SUi2),ns = 0] 
[50(3), ns = 0] 
[5f/(3),n5 = 3]. 



(3.3) 



This shows that there is also a compression within the sfermion mass spectrum, as the 
ratio of the squarks to the slepton masses is decreased in the extended models compared 
to the MSSM in the ttiq = limit, or more generally for any given value of niQ. This is 
because of the increased relative importance of the contribution to scalar masses from large 
renormalization scales where all of the gauge couplings and gaugino masses are larger. 

Despite these compressions in the gaugino and sfermion sectors considered separately, the 
combined sparticle spectrum in the extended models is stretched rather than compressed 
compared to the MSSM. Comparing eqs. (13.11) and (13.31) . one observes that in each of the 
extended models, the bino is much lighter than the lightest slepton, and so the lightest 
supersymmetric particle (LSP) will be a neutralino. This is in contrast to the well-known 
fact that the ttiq = scenario in the MSSM problematically predicts a stau as the LSP. The 
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supersymmetry-breaking flavor problem thus can be naturally solved by taking mp <^ 1^1/2 
in the extended models without running into the difficulties found in the gaugino mass 
dominated MSSM. 

Now consider the soft supersymmetry breaking masses for the new particles. Assuming 
gaugino mass unification, the Gx gaugino is heavier than all of the MSSM gauginos in 
the minimal Gx = SU{2) and SU{3) cases, but it is lighter than the MSSM gauginos if 
Gx = 5*0(3). In terms of the unified gaugino mass parameter mi/2, one finds at Q = 1 
TeV: 



Mj^ = mi/2 X < 



1.30 [SU{2), ns = 0], 

0.051 [50(3), ns = 0], (3-4) 
0.65 [SU{3), ns = 3]. 



[Compare eq. (13.11) .] The scalar members of the D, D and L, L multiplets get RG contribu- 
tions to their soft masses from both Gsm and Gx gaugino loops. Therefore, they are heavier 
than their MSSM counterparts with the same gauge quantum numbers. For the case where 
mi/2 dominates, one finds approximately for the soft masses ntf^ = and mj^ = mj., again 
at g = 1 TeV: 



) = mi/2 X < 



;i.97, 1.40) [SU{2), ns = 0], 

;i.22, 0.65) [50(3), ns = 0], (3-5) 

;i.68, 1.51) [5f/(3), ns = 3]. 



Also, for the minimal SU{3) case with ns = 3, one finds that nig = = 1.48mi/2- This 
is only slightly lower than m^, and m^, because most of the RG contribution to these soft 
masses comes from X loops in this case, which are the same for all of the new scalars. 

The qualitative features of the above results are illustrated in Figure [2j The soft masses 
for the gauginos, the first-family sfermions, and the new scalars are shown for Q = 1 
TeV. For purposes of comparison, mi/2 is chosen so that the heaviest MSSM squark, q, 
has the same mass in each of the four cases. The mass spectra in the extended models 
are readily distinguishable from the usual "mSUGRA" case parameterized by mi/2, "'tiq; ^o- 
This is because obtaining such a large ratio of scalar masses to gaugino masses in mSUGRA, 
would require a large mg, which in turn would lead to a much more compressed scalar mass 
spectrum. In contrast, the extended models are characterized by relatively heavy scalars 
which nevertheless maintain a significant hierarchy between squarks, left-handed sleptons, 
and right-handed sleptons, especially in the Gx = 5*0(3) and 5*^7(3) cases. 

In order to keep the discussion bounded, I will not give detailed results on the extended 
models with more singlets [ns > for Gx = SU (2) and 5*0(3), and ns > 3 for Gx = SU{3)]. 
However, the following qualitative features are notable. First, at one loop order, the presence 
of additional Gsm singlets does not affect the RG running of MSSM-field soft terms, so the 
effects are rather mild on the gluino, wino, bino, and MSSM squark and slepton masses. 
Second, increasing ns will decrease both the gx gauge coupling and the Gx gaugino soft 
mass at lower RG scales. Therefore, the Gx gaugino mass will be smaller compared to 
the MSSM gaugino masses Mi, M2, and M3 than in the cases shown in Figure O Also, the 
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FIG. 2: Comparison of the soft supersymmetry-breaking mass spectra following from a unified 
input gaugino mass that dominates at the unification scale, for the MSSM and the minimal 
models with Gx = SU (2), 50(3), and SU (3). The labels 1, 2, 3, X refer to running gauginos masses 
for U{1)y-, SU{2)i, SU{3)c and Gx, respectively. The labels q,u,d,i,e refer to MSSM first family 
squark and slepton soft masses, and the soft scalar masses for vectorlike supermultiplets are labeled 
with symbols D, D, L, L, and (for the Gx = SU{3) case with ns = 3), S, S. The value of is 
chosen so that the heaviest MSSM squark, q, has the same mass in each of the four cases. 



soft masses = m-^ and mj^ = rrij- will become relatively smaller, tending towards the 
MSSM squark and slepton masses and respectively. The soft masses for rrig = will 
also decrease for larger ns, although they are always heavier than M^, which decreases faster 
for larger ns- For ns > 0, the Yukawa couplings k, k' can also come into play, decreasing 
the L, L, S, and S soft supersymmetry breaking masses. 



IV. THE fi PARAMETER AND THE LITTLE HIERARCHY PROBLEM 

The supersymmetric little hierarchy problem is a subjective but inspirationally important 
puzzle which questions the naturalness of viable model parameters. The essence of it is that 
once one applies constraints from the non-observation of a Higgs boson and of superpartners 
at both LEP2 and the Tevatron, the actual value of might be considered surprisingly 
low for generic soft supersymmetry breaking parameters. 

The largest loop correction to the mass in the MSSM is given by (in the decoupling 
limit m\o > "^^o)- 



mlo = m| cos^(2/3) + ^ sin^/? y'f 
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2 \2 



(4.1) 



where m^^ and m^^ are the top squark masses, and q- and si are the cosine and sine of the 
top squark mixing angle 9i. Now in the models discussed in and this paper, adding in 
the effects of the Yukawa coupling k one finds the further estimated correction in the case 
/Wl ~ /is ~ Mp with heavier scalars with masses of order Ms 



Amlo 



N 

47]-2 



sm 



73 



1 
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12M| 



(4.2) 



Here x = Mg/Mp is the ratio of the average new scalar and new fermion masses in the 
L, L, S, S sector, and Xk = ak/k — fi cot /3 is a mixing parameter for the scalars. The largest 
possible contributions come from the maximal (fixed-point) values of eq. (12. 25 p . As was 
pointed out in ref. jsl, for Gx = SU{3) eq. (14. 2 p is enough to raise the Higgs mass by tens 
of GeV, depending on the details of the fermion and scalar masses in the new sector. 

From the point of view of the supersymmetric little hierarchy problem, even raising 
the Higgs mass by a few GeV is potentially helpful. However, one must also consider the 
effect of the new sector on the scalar potential. The minimization of the Higgs potential in 
supersymmetry results in: 
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2 I 2 

+ m 



]__d_ 

iL, dv,. 



Ay + C(l/tan2/3), 



(4.3) 



where AV is the radiative part of the effective potential, with tan/3 = Vu/vd = {H^) 
and Vu treated as a real variable in the partial differentiation. In general, without further 
theoretical structure, ^ and have no reason to be related, since /x is a supersymmetry- 
preserving parameter and Tn^^ is supersymmetry-breaking. In the MSSM with generic 
parameters, one finds that —m'jj^ tends to be much larger than m|, and eq. (14. 3 p seems to 
imply a percent-level fine-tuning of the difference between and rnj^^. 

It is not possible to rigorously quantify fine tuning, since there can be no such thing as 
an objective measure on parameter space. Nevertheless, qualitative trends can be identified. 



and an obvious approach is to consider models with smaller predicted values of —tn'jj^ at the 



nin^, because then the fractional 



weak scale to be more likely than those with very large 
tuning required between it and will be less. This in turn means that smaller values of 
\n\ are more likely than very large values, since this is determined by eq. (14. 3p . 

With this in mind, it is interesting to consider how the MSSM and its extensions, and 
variations of the most popular models of supersymmetry breaking, affect the weak-scale 
predictions for —m^j^. For example, in the MSSM with tan/3 = 10 and rrit = 173.3 GeV, 
one finds from RG running at Q = 1 TeV in terms of the GUT-scale input parameters mi/2, 
Aq and mg: 



= 1.65m2/2 - 0.40mi/2Ao 



0.11^2 



0.022m^. 



(4.4) 



This formula shows that — m|^^, and therefore and therefore the level of fine-tuning 
required, increase with the gaugino squared masses. In extended models the gaugino masses 
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at the unification scale have a varying relationship with the gaugino masses at the weak 
scale, which are more closely related to the physical masses, so it is useful to reformulate 
this in terms of the running gluino mass parameter M3 also evaluated at Q = 1 TeV: 

- mjj^ = 0.32M| - O.I8M3A0 + 0.11 Al - 0.022m^, [MSSM]. (4.5) 



In fact, most of the dependence on the gaugino masses comes from the gluino mass [39|, 
so this formula is approximately valid even for moderate deviations from gaugino mass 
universality. One can note that for a gluino mass of order 500 GeV, and small |^o|) is 
only of order (280 GeV)^, so that the tuning needed to get m| in eq. (14. 3 p is of order 5%. The 
problem is that (although there is still considerable variation among models, particularly 
for large \Ao\) lower values of M3 typically give a prediction for m/jo that is smaller than 
114 GeV, and higher values of M3 require even more delicate cancellation between —m'jj^ 



and . The "focus point" region |40|, |41| occurs due to the small negative coefficient of 
ttIq in eq. (14. 5p . which allows a cancellation between the M| and ml terms for very large 
ttIq, leading to a small value of —mfj^ and therefore small However, this also can be 
judged to be fine-tuned, as the large value of rriQ has to be finely adjusted, given a value of 
the ostensibly independent parameter M3. 

We can now compare with the situation for the models in the present paper. In the 
minimal SU{2)x model with ns = 0, one finds instead from the RG running: 

- m^^ = 1.07M| - O.44M3A0 + 0.12AI - 0.16ml, [SU{2)x, ns = 0]. (4.6) 

again for —m\j^ and M3 evaluated at Q = 1 TeV, assuming gaugino and scalar mass univer- 
sality, tan/3 = 10 and = 173.3 GeV. (The coefficients change, but not very radically, for 
larger ns-) The larger coefficient of M| indicates that this is naively even more fine-tuned 
than the MSSM. However, this effect is not without compensation; as Figure [2] shows, one 
does not need as large a gluino mass to get large squark masses, which in turn lead to 
large positive contributions to m^o from eq. (14. ip . Also, the larger negative coefficient of m^ 
means that the analog of the MSSM focus point region occurs at much smaller values of m^ 
in this extended model. Since there is no such thing as an objective quantitative measure 
of fine-tuning, I choose not to attempt to make a definitive statement beyond observing the 
competing factors just mentioned. 

For SO{?>)x with ns = 0, the analogous formula becomes: 

- mjj^ = 3.36M^ - 0.91M3ylo + O.lSAl - 0.26ml, [S0{3)x, ns = 0]. (4.7) 

Similarly, for the minimal viable SU{3)x model with ns = 3 and small Yukawa couplings 
k, k' ~ 0, the analogous formula becomes: 

- m]j^ = 5.97M| - I.28M3A0 + O.UAl - 0.28ml, [SU{3)x, ns = 3]. (4.8) 

In both of these cases, the situation again seems subjectively worse with respect to fine- 
tuning than the MSSM, due to the much larger coefficient of M|. As shown in Figure |2l one 
does naturally get much larger squark masses for a given M3, again leading to larger radiative 
corrections to mf^o. However, with these large coefficients to M3, the direct (but model 
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dependent) constraints on the gluino mass from Tevatron come into play in a significant way. 
Again, there is a chance for more cancellation between the gaugino and scalar contributions 
due to the negative coefficient of mg. 

It is also interesting to consider the situation for > 1 with large Yukawa couplings 
near the fixed points of eq. f l2.25p . For SU{2)x with ns = 1 and the fixed-point value 
k = Ajfixcd = 0.88 at Q = 1 TeV, one finds: 



m 



\^ = 1.28M| - O.I5M3A0 + 0.04^2 + 0.72ml [SUi2)x, ris = I, k = kf,,,^].{A.9) 



Here the coefficient of the gaugino mass squared is even larger than for k = [compare 
eq. (14. 6p ]. and the coefficient of the scalar squared mass rriQ is large and positive, eliminating 
the possibility of cancellation to achieve a smaller —m\j^. Similar results obtain for SO{?))x 
with ns = 1 and k = fcfixed = 0.76 at Q = 1 TeV: 

- = 3.46M| - O.27M3A0 + 0.05AI + 0.89ml, [S0{3)x, ns = I, k = Affixed], (4.10) 
and for SU{3)x with ns = 3 and one k = fcfixod = 1-32 at Q = 1 TeV: 

- mjj^ = n.OMi - 0.21M^Aq + 0.02^1^ + O.Tlm^, [^f/(3)x, ns = 3,k = A;fixed]. (4.11) 

Therefore, even though the fixed-point Yukawa coupling can give large positive contributions 
to 771^0, there is a quite detrimental effect on the fine-tuning needed to obtain the observed 
mz in models that have heavy enough gluinos (and charginos) to have evaded discovery at 
the Tevatron and LEP2. Similar effects have been noted before in the case of vector- like 
fermions without an additional gauge group in refs. [lo[, ll||, (l3| . 



Qualitatively, the model with SU{2)x and no new Yukawa coupling seems to be the least 
fine-tuned of the extended models. Adding new Yukawa couplings, despite increasing mhO, 
does not clearly alleviate the little hierarchy problem, and arguably make it much worse, 
especially in the cases of S0{3)x and SU{3)x- 



V. COLLIDER PHENOMENOLOGY OF THE QUIRKS 

In this section, I will consider some features of the phenomenology of the quirks in the 
models discussed above, following for the most part general ideas and results from refs. 18 



19[, [20|, and [21[. For simplicity, I will consider only the fermions from the D,D, L,L, 
multiplets, and not their scalar partners. This is because supersymmetry breaking effects 
provide for the scalars ( "squirks" ) and the Gx gaugino to have much larger masses, making 
them less immediately relevant for collider searches. (Even if they had the same masses, 
squirks would have much smaller production cross-sections than fermionic quirks. The Gx 
gauginos will not be produced directly in tree-level processes at colliders at all.) When 
produced, the squirks will decay promptly to quirks and MSSM gauginos. The Gx gaugino 
can undergo a three-body decay to a quirk, antiquirk and MSSM gaugino, if kinematically 
allowed. In this section, I will use the same symbols for the fermions as for the chiral 
supermultiplets to which they belong. 
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I will also assume for simplicity, and motivated by the results of the previous section, that 
the mixing of the Gsm singlets S, S with the doublets L, L due to the Yukawa couplings k, k' is 
small in most of the following. This implies that 5*, S decouple from collider phenomenology. 
The charged {L~,L^) and neutral {L^,L ) fermions form two Dirac fermion-antifermion 
pairs, each with tree- level mass fii- However, radiative corrections split the masses slightly, 
with Am = rriL- — rrij^o always positive. (If present, the Yukawa interactions k, k' that cause 
mixing with S, S would increase this splitting, so the lightest non-colored fermion is always 
neutral.) One finds Am > 270 MeV for fii > 100 GeV, with Am approaching 355 MeV 



asymptotically for large fii 42|. This means that the decays 



^ T^i+ue, (5.1) 
L- L'^rVi, L°7r-, (5.2) 

(and decays to more pions or other SM hadrons if non-zero k, k' increase Am) mediated by 
the W boson are always kinematically allowed, and will occur with decay lengths of order 



centimeters [42] due to the small available kinematic phase space. Therefore, the lifetime of 
L^,L~ is large compared to other processes to be discussed below; in particular they will 
form quirk-antiquirk bound states and annihilate before they decay. In the simplest scenario 
(barring additional couplings to be described in the next paragraph), the neutral quirk Dirac 
fermions L°,L are completely stable, as are the colored quirks D,D with charges ±1/3. 
Note that none of the quirks can mix with the Standard Model fermions because of Gx 
conservation, so the lightest quirk is always stable. Such stable fermions could present a 
challenge for the standard cosmology with a high reheat temperature but need not be a 



disaster 43 , 44 , 18 



If ris' > 1 and the pairs L, L and S, S have the opposite matter parity from each other, 
then the Yukawa couplings k, k' are forbidden, but the superpotential term 

W = XeSle (5.3) 

is allowed, with i an MSSM SU{2)l doublet lepton. This provides additional possible decay 
modes L — )■ S£^ and L" — )■ Si~[, if \ fill > \f^s \ +'^/ so that these decays are kinematically 
allowed, or alternatively with the sleptons off-shell. These decays are not automatically 
kinematically suppressed, and so could happen promptly before quirk-antiquirk annihilation 
occurs. The fermions D and D may also have additional decay modes, if the possible 
superpotential terms 

W = \LDq + \-^Dd, (5.4) 

are present. The first of these terms is only allowed if Gx is either SU{2) or S0{3), and 
if L,L have the opposite matter parity of D,D (if matter parity is conserved). It permits 
squark exchange to mediate the decays 

D L^dL, L-ul, (5.5) 

D l?dl, T^ul, (5.6) 
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with the MSSM squarks possibly off-shell due to kinematics. The second term in eq. (15. 4p 
is only allowed if > 1 and 5", 5" have the opposite matter parity of D. Then MSSM 
squark exchange can mediate the decays 

D SdR, D Sd*ji, (5.7) 

again with the squarks possibly off-shell. Whether these decays can be important depends 
on the kinematics as well as the size of A^, Xg, and A^. For simplicity, they will be assumed 
to be absent or at least too small to make a difference below, except where noted otherwise. 
The important direct pair-production processes for the new fermions are 

pp DD, (5.8) 

pp Z^*\-f^*^ L^L~, 1^L'>, (5.9) 

pp W+^*^ T^L^, (5.10) 

pp W^-W ^ X°L-, (5.11) 

for the LHC, with the obvious substitution of pp for the Tevatron. Pair-produced quirks with 
masses much larger than A will move apart from each other with typically semi-relativistic 



speeds, and as described in [l8|, will be connected by Gx flux strings with tension a. From 



the lattice, there is an estimate (see Table 7 and eq. (11) of |45j ) : 



-a^K^A''' ^''<^>- (5,12) 
I 1.86 5f/(3)x, 

so that the maximum string length in a given hard scattering event with kinetic energy Ai? 
in the center-of-momentum frame is 

A^ f AE \ fkeV\^ 
L = ^ 6 mm — — . 5.13 

VlOOGeV; ' 

Therefore, the lengths of such strings, although much larger than A~^, will typically be less 
than 1 mm for A greater than a few keV. From Table [H one finds that the quirky flux strings 
will be microscopic for SU{2)x with ns <3 and for SU{3)x with ns < 9, assuming that gx 
is unified with the SM gauge couplings and all singlets charged under SU{N)x are heavier 
than A. 

For the case of Gx = 5*0(3), the situation is quite different, because from Table [T] the 
confinement distance scale A~^ is literally astronomical, at least of order the Earth's orbit 
around the Sun even in the minimal model. The quirks in this case are essentially free 
particles with multiplicity 3 times larger than expected from their SM quantum numbers. 
Note that even if one rejected the unification of gx with the SM gauge couplings in this 
model to arrive at a much larger A, the fact that the supermultiplets are in the adjoint 
representation of the Lie algebra means that they would not form stable flux tubes of the 



type discussed in 18|] when pair-produced, even if any S, S fields are heavier than A (so 



n = 0). Instead, pair-produced particles charged under S0{3)x would each bind to a gauge 
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boson to form two stable Gx-singlet states with size of order A~^, allowing the flux tube 
to break. Although the new fermions behave like free stable particles when pair-produced 
at colliders, the fact that they will come in three-fold exactly degenerate multiplets will in 
principle allow a determination of their nature from their production cross-sections. In the 
simplest case, D, D and L , L° will be absolutely stable, with L^, L~ having decays to L ,LP 
via soft pion or lepton emission as discussed above. The L , are only weakly interacting 
and thus invisible, but known collider search strategies 46| for stable strongly interacting 
particles apply for D,D. However, as noted above, D, D may be able to promptly decay 
according to eqs. (]5.5p -( l5l7|) . depending on both kinematics and the allowed superpotential 
terms. If so, then the signatures will always contain E^^^^, and will resemble those for 
ordinary MSSM squarks. 

For the remainder of this section, consider the cases of SU{2)x and SU{3)x, with the 
Gx confinement scale less than the masses of the quirks that have SM gauge interactions, 
and stable microscopic flux strings joining the quirk-antiquirk pairs. The quirk- ant iquirk 
pair will then form an exotic bound state with invariant mass given approximately by the 
total center-of-momentum energy of the hard partonic scattering that produced them. This 
quirk-antiquirk string state can lose energy either by G^-glueball emission, by radiation of 
many soft photons, or in the case of the DD state by radiation of numerous soft pions, a 



"hadronic fireball" 18|. The large multiplicity of soft pions or photons may be detectable 
as anomalous "underlying events" 18|, |21[ that accompany the hard scattering production, 
and may be used as an additional tag to dramatically reduce backgrounds. 

If the quirk and antiquirk lose most of their initial relative kinetic energy before annihi- 
lating, they will briefly form a "quirkonium" bound state which then decay^ to two or three 
hard partons with invariant mass peaked at twice the mass of the quirk Jl8|. Alternatively, 
however, the neutral and colorless quirk and antiquirk states might [l8| have a prompt 
annihilation before they can lose enough energy to form a low-lying quirkonium state. In 
that case, the final states will have a broad distribution of annihilation products, which will 
therefore be much harder to discern above hadron collider backgrounds. It is difficult to 
estimate in advance what proportion of the events will fall into these two categories, due to 
the non-perturbative nature of the energy loss mechanisms, which do not have direct analogs 
in experimentally known hadronic physics. 

For the weakly interacting quirks, and for the strongly interacting quirks if A > Aqcd? one 
might suspect the non-perturbative interactions by which the quirk-antiquirk string state 
loses energy to be dominated by Gx-glueball emission. However, this is quite uncertain, and 
can be suppressed or even eliminated by kinematics if A S> A pr;]-) . The masses of the Gx 



glueballs have been estimated by lattice computations [47|-[49| 
lightest two glueball states with J^*^ = 0"'""'' and 2"'""'": 



45 , with the results for the 



mo++ 
m2++ 



Q.7A- 
9.6A 



MS' 



MS- 



(5.14) 
(5.15) 



There are other heavier glueball states 0"'""'"*, 
only there are also states with odd G, , 



3+-, 



2-+, 

0+- 



0-+*, 
1", 



2-+* 
2— 



, and, for SU{3)x 
3 , with masses 



ranging up to about 3mo++. As can be seen from Tables HI and HTl in the case of SU{3)x 
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these glueballs should have masses in the hundreds of GeV range for the minimal case of 
ns = 3 and so could be comparable in mass or even heavier than the lighter quirks, and 
should be in the tens of GeV range for ns = 4. This would prohibit energy loss of the quirk- 
antiquirk flux string states into Gx-glueballs. For the other cases listed in Table HI decays 
of the flux strings to G^-glueballs should be allowed, but perhaps kinematically suppressed, 
leading to considerable uncertainty in the number of G^-glueball states emitted and the 
likelihood of the quirk-antiquirk string state to lose most of its energy before annihilating. 
It is also possible that a few Gx-glueballs will be produced in the original hard scattering 
production. 

If produced, the detection of Gx-glueballs is problematic. Their decay widths can be 



estimated for SU{3)x using eqs. (17), (23) and (30) of ref. [231] (see also ref. |26|]) with 
matrix elements from eqs. (38) and (62) of ref. (5o| : 

r(0++ ^ gg) = 360 4a7/o, (5.16) 

T{2++^gg) = 0.12 «|A7/i8,, (5.17) 

TiO-^^gg) = 24a|AV/o. (5.18) 

The results for SU{2)x should be comparable and slightly smaller. This leads to proper 
decay lengths for Gx-glueballs of order 



/0.2\V /^D ^^fGeVV J 0.14 meters (for 0++), 
^"=UJ (i^^) V^J "1 400 meters (for 2++), (5.19) 

2 meters (for "'"). 

If A ^ 1 GeV as expected for the minimal SU{2)x model, a sizable fraction of the O^"*" 
decays might occur within the detector, but only if /i^) is less than roughly 150 GeV, which 
may be ruled out already by Tevatron data (see below). For larger fin or smaller A, the 
decays of the Gx-glueball will occur outside of the detector and will be invisible. For much 
larger A as occurs in the SU{3)x model with 77,5 = 3 or 4, the decays may occur within 
the detector for any /i^j, but then the production of Gx glueballs in the flux-tube energy 
loss processes is likely irrelevant anyway due to kinematic suppression or prohibition. Even 
if the Gx glueballs are produced and decay promptly, the main decay is likely to a pair of 
gluons, and the resulting dijet mass peak signal from these decays will have to compete with 
a huge background from QCD. To have a significant branching fraction to 77, which has 



much smaller backgrounds, one can take fiL < fJ'D, with a leading-order estimate (23|, [26 



BR(0++^77) _8a' fl Sfi^V 



For example, with ~ 0.2, this ratio is of order 0.001 for /iL = /^d? but it rises to about 0.18 
if yUi/ fiD = 0.5, and is greater than 1 if fii/f^D is less than 0.4. The non-perturbative nature 
of the Gx glueball production mechanisms means that the diphoton signal strengths, if any. 
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are extremely difficult to estimate even roughly, but to have even a hope of observation would 
seem to require Hl < IJ'D and A of order a few GeV (not too small for cr to be large, but not 
too large for Gx-glueball production to be kinematically suppressed). Nevertheless, given 
the uncertainties involved, this possibility highlights the general importance of searching for 
narrow diphoton peaks at large invariant masses at the LHC; this type of signal could also 
arise not only for the classic diphoton signal for a low-mass Higgs scalar boson, but also for 
stoponium 51, 52] or for Kaluza-Klein gravitons in theories with low-scale gravity [53|. 

Probably the most optimistic scenario for detecting the quirks occurs in the case that DD 
are strongly produced at a hadron collider and manage to lose most of their initial relative 
kinetic energy stored in the Gx flux tube by radiating soft pions and/or Gx glueballs, 
arriving at a low- lying quirkonium state with mass ~ 2/i/) before finally annihilating in a 
color-singlet S-wave '^^^^Lj = ^Sq (r/) or ^Si (ip) state. The most promising channel for 
detecting the quirkonium peak is a'^ti~. The relevant annihilation decay widths for a ^Si 



state can be inferred from refs. 54 and l2d 



'l + 2Rf)(ef- 



9) 



Cvi/(1 — Rz] 



a^ell3lyl + 20Rw + l2Rl,. 



(1 - Rz)' 



T{^p-^ggg) = ^in'-9)To 



7V3 



9)ro 



To, (5.21) 

(5.22) 

(5.23) 
(5.24) 



where e/ = (2/3,-1/3,-1,0) and A^/ = (3,3,1,1) and gf = (1/4,-1/4,-1/4,1/4) for 
/ = {u,d,e,h') respectively, and gj = gf — e/s^, and sw and cw are the sine and cosine 
of the weak mixing angle, and eo = —1/3, and Ri = mf/M'^, where M ^ 2/i/) is the 
quirkonium mass, and /3j = y^l — ARf , and Fq is a common normalization proportional to 
the square of the wavefunction at the origin. The Gx gluon is represented by X. Note that 
final states gg, ZZ, Zj, 77 and XX do not occur in ip decays. The final states Zgg, 'jXX, 
ZXX do occur, but with branching ratios that turn out to be very small. For N = 2, the 
decay to three SU (2) y gauge bosons vanishes due to the A^^ —4 factor, and for SU (3)x with 
A^ = 3 the XXX decay will be kinematically forbidden or at least highly suppressed by the 
large masses of the Gx glueballs that would have to be the final result of Gx-hadronization. 
Therefore final states involving Gx glueballs should not play a significant role in quirkonium 
decays. 

For ^5*0 states, the dominant decay is to gg or XX, and // does not occur at all at leading 
order. If we assume that the spin state is randomized by the non-perturbative processes that 
lose the initial relative kinetic energy, so that ip and 77 states are populated in the ratio of 3 to 
1, then the branching ratio of quirkonium to leptons should be given by 3/4 of the branching 
ratio indicated by eqs. ( KT[^ -( E?2^ . Numerically this yields^ BR{DD j^^l^') = 0.093 



^ The estimate in section 5.6 of 



18] is parametrically different, and numerically smaller by a factor ~5. 
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FIG. 3: The production cross-section for DD 
in pp collisions at t/s = 14, 12, 10, and 7 
TeV, and in pp collisions at -y/i = 1.96 TeV, 
obtained at leading order using CTEQ5L0 



57l | parton distribution functions (PDFs) with 
Q = fJ-D- 



very nearly independent of the mass. It should be noted that this branching ratio does 
not apply to prompt annihilation of the quirks before they have settled into a color- and 
Gx-singlet quirkonium state; that branching ratio will be much smaller, and will not lead 
to a sharp dimuon peak, and so leads to a more pessimistic case. 

In the most optimistic case that most of the DD states annihilate after losing most of 
their excess energy, there are good prospects for detection at hadron colliders, because the 
signal production is strong and peaked in invariant mass, while the dominant background 
is electroweak (Drell-Yan) and diffuse. The total production cross section at the Tevatron 
and at various LHC energies is shown in Figure |3l The CDF collaboration has published 
55[ a limit on cross-section times branching ratio for new states that decay to /x'''/^", based 
on 2.3 fb~^ of pp collisions at the Tevatron. Comparing the relevant spin-0 limit from 
Figure 3 in 55| to the results shown in Figure [3] of the present paper and using the estimate 
BR(yU"^/i~) = 0.093 from above, I obtain the lower mass bound fio > 375 GeV in this 
optimistic case. 

At the LHC, the invariant mass resolution for high-mass dimuons should be of the order 



of 5% (56| for the CMS detector. Therefore, as a rough estimate of the discovery reach, I 
consider a mass window from 0.9M to I.IM where M ^ 2/xd is the quirkonium mass, and 
require that S/ y/B exceeds 5 in that window, where S is the number of signal events (which 
is also required to exceed 10) and B is the expected number of Drell-Yan background events. 
The Drell-Yan background cross-section is shown in Figure HI Trigger and detector efficien- 
cies are not included, but these are expected to be very high for high-mass dimuon events, 
and the QCD i^-factor for the signal is not included. Dimuon backgrounds from sources 
other than Drell-Yan can be suppressed by requiring no extra hard jets or missing energy. In 
the following, I will again assume a spin-averaged BR(/i"'"/i~) = 0.093 for the signal. There 
is also a potential confirming signal from annihilation to e~^e~ , with an invariant mass peak 
that is similar but wider and smaller due to detector resolution and efficiency effects. 

For a 1 fb~^ LHC run at ^/s = 7 TeV, the signal cross-section in Figure |3] yields 20 
expected dimuon events for hd = 500 GeV, and as shown in Figure H] there is about 1 
background event expected in the corresponding mass window M{fi'^fi~) = 1000 ±100 GeV. 
Requiring 10 signal events, the discovery reach is up to about fio = 550 GeV. 
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FIG. 4: The leading-order differential produc- 
tion cross-section for the background, 
through 7* and Z*-mediated (Drell-Yan) pro- 
cesses in pp collisions at ^/s = 14, 12, 10, and 7 
TeV, and in pp collisions at ^/s = 1.96 TeV, as 
a function of the invariant mass of the 
pair. CTEQ5L0 PDFs were used with 



For LHC pp collisions at a/s = 14 TeV, the signal cross-section times dimuon branching 
ratio for fio = 800 GeV is 15 fb, with a background level in the mass window M{fi'^fi~) = 
1600 ± 160 GeV of 0.8 fb. Therefore, discovery may be possible in this case with 1 fb~^. The 
mass reach is essentially determined by the number of signal events, since the background 
levels in the high-mass windows are small. In the same way, with 10 fb~^, I estimate the 
10-event discovery reach to he up to fio = 1200 GeV, and for 100 fb~^ up to about fj,D = 
1600 GeV. 

In a more pessimistic scenario, the quirk and antiquirk may usually annihilate before they 
can settle into a low-lying color-singlet quirkonium state. The branching ratio to dileptons 
will be severely reduced in that case because there are color octet as well as color singlet 
decay states available, and the remaining dimuons will be distributed over larger invariant 
masses. If one supposes that only 10% of the DD pairs that are produced will settle into a 
low-lying color-singlet quirkonium state before annihilation, and uses only the dimuon events 
from this quirkonium peak, then the signal cross-section before BR(yU^yU~) is effectively ten 
times smaller than shown in Figure [3l The limit from comparing to the CDF bound on 
cross-section times branching ratio (Figure 3 in 55|, based on 2.3 fb~^) results in /x^ > 180 
GeV. I estimate that the expected reach from a 1 fb~^ LHC run at ^/s = 7 TeV in this more 
pessimistic case is roughly fi^ = 350 GeV, for which about 17 dimuon signal events and 6 
background events would be expected in a mass window M(/i+/i~) = 700 ± 70 GeV. For 
LHC runs at ^/s = 14 TeV with (1, 10, 100) fb~^, I similarly estimate that the discovery 
reach for DD that annihilate at least 10% of the time from color-singlet 5'-wave quirkonium 
would extend to about hd = (500,800, 1100) GeV. 

In the case of the non-colored quirks L^, L", L ,L^, the production rates are electroweak, 
and the energy loss rate for the quirk-antiquirk bound by the flux string is much lower 



18[. The most promising signal may come from the production of the quirk-antiquirk 
states with a net ±1 charge, as in eqs. f l5.10p and (15. lip , because charge conservation then 
prohibits the subsequent prompt annihilation to invisible Gx glueballs that may occur in 
the case of neutral bound states. The analogous case for fractionally charged squirks in 



"folded supersymmetry" was proposed and studied in 19|. The excess energy from the hard 



production will be radiated away in the form of Gx glueballs or soft photons, hopefully 
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FIG. 5: Branching ratios for charged quirkonium (a low- lying or L bound state) in a 

(left panel) or a ^Si (right panel) state. 



allowing the quirk and antiquirk to finally annihilate when nearly at rest in a charged 
quirkonium bound state. The annihilation is strongest in an S'-wave state. I will again 
assume that spins are randomized by the energy loss process, so that '^Si and states are 
populated in the ratio of 3 to 1. The branching ratios for such states have been computed in 



ref. [20|, and are shown in Figure [5]for the present case of constituent quirks with charges ±1 
and 0. The ^5*0 state decays predominantly into W'-f, with an invariant mass of nearly 2/1^, 
and therefore a hard photon. (A somewhat smaller branching ratio to W'^ was obtained in 



ref. 20[ for a case with fractionally charged constituent quirks.) This state may therefore 
be searched for in the £^7 + E^^^^ channel at hadron colliders, as suggested in the similar 



squirk case of ref. 19 



The combined''- production cross-sections at the Tevatron and at various possible LHC 
energies for the charged quirk-antiquirk combination are shown in Figure [6l These cross- 
sections are about an order of magnitude larger than for fractionally- charged scalar quirks 
(as studied in ref. [lij) of the same mass. Partly counteracting this, one might expect that 
only about 1/4 of fermionic quirk-antiquirk production will end up in a ^Sq state that can 
annihilate to W^7, rather than a ^Si state that decays mostly to jets or a single lepton plus 
neutrino. Thus, the effective branching ratio of charged quirkonium should be about 0.2 for 



W^, a factor of 3-4 smaller than used in 19j. The net effect is that the total production 



cross-section times branching ratio for 1^7 should be a factor of 2-3 times larger, for a given 



quirkonium mass, than in the study of ref. 19 



The largest background is from Standard Model W'-f production, which features a rapidly 
falling tail at high photon px- In contrast, the signal from ^Sq quirkonium decaying to W'^ 
should have a photon px distribution that is approximately flat, with an endpoint near 

— m'^/4fiL in the idealized case that the transverse kick to the quirkonium is small. The 



^ The charge +1 combination is produced more often than the charge —1 one at the LHC, as usual. 
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FIG. 6: Combined leading-order production 
cross-section for and L L~ in pp colli- 

sions at the LHC for ^/s = 14, 12, 10, and 7 
TeV, and in pp collisions at the Tevatron with 
= 1-96 TeV. CTEQ5L0 PDFs were 
used with Q = fJ-L- 



relevant photon pj- distribution has been studied at Tevatron by both CDF 58|] and D0 (59 



where it was found that the data is described well by the SM W'-f and other subdominant 
backgrounds including Wj with the jet faking a photon and Z7 with one lepton from the Z 
missed. At hadron colliders, a W'j mass peak can in principle be reconstructed if one assumes 
that the observed i?™'^^ in the event is due to the neutrino from the leptonic W decay, but 
this is subject to the considerable uncertainty in how much missing energy is actually due 
to missing Gx-glueballs radiated from the initial state, as well as from the underlying event, 
additional jets, or from mismeasurement. However, the discovery potential may be greatly 
enhanced because one can also look for a large number of soft photons radiated as the quirk- 
antiquirk flux string loses energy, forming an anomalous "underlying event" with distinctive 
character. The resulting complications are beyond the scope of the present paper, but have 



been discussed in the analogous case of fractionally charged colorless squirks in [19|, Ul 
The search for W'-f candidates a with large photon pt and a possible peak in invariant 
mass or transverse mass, in combination with an anomalous underlying event used as a 
background-reducing tag, may well be the best hope to detect the quirks in these models. 



VI. CONCLUSIONS 

Extensions of minimal supersymmetry with an extra non-Abelian gauge group and quirk 
supermultiplets maintain two of the hallmark successes of the MSSM: compatibility with 
perturbative gauge coupling unification and with constraints on precision electroweak ob- 
servables. Natural mechanisms can put the quirk fermion masses at the TeV scale or below. 
It follows that requiring the unified gauge couplings to be perturbative, so that low-scale pre- 
dictivity is not lost and the apparent unification of gauge couplings is not just an accident, 
the gauge group under which the new vector-like particles transform in the fundamental 
representation must be either SU{2), SU{3), or 5*0 (3). The presence of a new non-Abelian 
gauge group has a dramatic effect on the superpartner mass spectrum, and allows the soft su- 
persymmetry breaking parameters to be dominated by the gaugino masses at the unification 
scale, while still having a neutralino LSP. 
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In the 5*0(3) case, the confinement length scale is so large as to make the new particles 
essentially free in collider experiments. In contrast, in the minimal versions of the SU(2) and 
SU (3) cases, the quirk- ant iquirk bound states produced at colliders will be microscopic. If a 
significant fraction of the colored quirk-antiquirk pairs produced at hadron colliders will lose 
most of their excess energy before annihilating as quirkonium, then there is significant reach 
in the dilepton mass peak channel. Even if this fraction is only 10%, Tevatron data that has 
aheady been analyzed should allow a limit of 180 GeV for the quirk mass to be set from the 
search for a resonance. In the optimistic idealized case that all of the quirk-antiquirk 

pairs annihilate from quirkonium, this limit should be about 375 GeV. In a 1 fb~^ LHC run 
at y/s = 7 TeV, the discovery reach could be as high as 550 GeV, and should extend well 
above 1 TeV with 100 fb"^ at y/s = 14 TeV. The color singlet quirks in these models can 
also be searched for as quirkonium W'-f resonances. In all cases, the quirkonium peak can be 
accompanied by an anomalous underlying event consisting of many soft pions or photons, 



which can significantly aid in making a discovery |18l . |21| . If low-energy supersymmetry is 
realized in nature, then it will be important to test the possibility that it is not minimal by 
searching for these events. 
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